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We consider entanglement for quantum states defined in vector spaces over 
the real numbers. Such real entanglement is different from entanglement in 
standard quantum mechanics over the complex numbers. The differences 
provide insight into the nature of entanglement in standard quantum theory. 
^ ' Wootters [Phys. Rev. Lett. 80, 2245 (1998)] has given an explicit formula for 

f^ ' the entanglement of formation of two qubits in terms of what he calls the 

f— ^ . concurrence of the joint density operator. We give a contrasting formula for 

On , the entanglement of formation of an arbitrary state of two "rebits," a rebit 

^^ I being a system whose Hilbert space is a 2-dimensional real vector space. 

o' 
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I. INTRODUCTION 



H ! One of the key distinguishing features of quantum mechanics, not found in classical 

^ I physics, is the possibility of entanglement between subsystems. The significance of this 

^1 phenomenon is now unquestioned, as it lies at the core of several of the most important 

achievements of quantum information science |jl[, such as quantum teleportation and 

^ . quantum error correction [0,§]. Yet can we say that we understand the distinction between 

^ ! those physical theories with entanglement and those without? It is difficult to claim such 

understanding, as our most well studied foil theory to date — namely, classical physics — is 

completely devoid of the phenomenon. This paper, in a small way, contributes to filling that 
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gap by studying entanglement in real vector spaces, where though there is entanglement, it 
is different from the entanglement in the complex vector spaces of standard quantum theory. 

Our specific objective here is to analyze entanglement for states of two "rebits," a rebit 
being a two-state system whose Hilbert space is defined over the field of real numbers. 
In particular, we give an explicit formula for the entanglement of formation of an arbitrary 
(mixed) state of two rebits. This formula is similar in some structural ways, but not identical 
to Wootters's formula for the entanglement of formation of two qubits [^]. The reason for our 
considering real vector spaces is just the reason described above: they provide an easy, well 
defined foil with which to compare standard quantum theory; indeed, they have already been 
fruitful in that regard 0. A quantum theory that uses real vector spaces is similar to, but 
not identical to the standard theory [^H- In general, we follow the philosophy of Weinberg 
in Ref. in the hope that this exercise will help identify those aspects of entanglement that 
are unique, those that are accidental, and those that are necessary to the standard theory. 

To state our main result, we must build some concepts and notation. Consider a bipartite 
composite system, made up of subsystems A and B. A density operator p^^ of the composite 
system, pure or mixed, is said to be separable if it can be thought of as arising from an 
ensemble of product states, i.e.. 



p"'=Ew"®pr (1-1) 



A separable pure state is itself a product state. The reason this definition is interesting is 
because a separable state can be created by procedures that are local to each subsystem, 
whereas a nonseparable state cannot be created by any local means. 

Taking the matrix transpose of any density operator relative to some orthonormal basis — 
this is the same as taking the complex conjugate in that basis — yields another density 
operator, i.e., another positive semi-definite operator with unit trace. Similarly, if a state 
of a bipartite system is separable, taking the partial transpose on system B in any basis 
also yields another density operator. If, however, taking the partial transpose leads to 
an operator that is not positive semi-definite, one can be sure that the original state was 



an entangled state. This is the partial transpose condition of Peres ||T0|. Unfortunately, 
for subsystems A and B of arbitrary Hilbert-space dimensions, the converse of the Peres 
condition is not true — the partial transpose of an entangled state can give another positive 
semi-definite operator. Thus the Peres condition does not provide a general criterion for 
testing entanglement. For 2x2 systems (two qubits) or 2 x 3 systems (a qubit and a qutrit), 
however, the Peres condition does provide a criterion for entanglement: a state of such a 
composite system is entangled if and only if its partial transpose is a nonpositive operator, 
i.e., has at least one negative eigenvalue |[TT[] . 

The chief resource-based measures of entanglement are the entanglement of formation 
and the distillable entanglement |T2| , p!3| . For two d-dimensional systems, the pure state 



I^) = ^EIS")®K")' (1-2) 

where \ej) and |ej) are orthonormal bases for the two subsystems, is maximally entangled 
in the sense that it can be used to teleport the state of another d-dimensional system. The 
degree of entanglement of such a maximally entangled state is logj d, the marginal entropy of 



each subsystem. Suppose that starting with m such maximally entangled states, one has a 
procedure, involving only local operations and classical communication between subsystems, 
for creating n copies of an arbitrary state p^'^ . The entanglement of formation, E{p^'^), is 
defined to be \0g2d times the asymptotic ratio m/n for an optimal procedure, i.e., one that 
has the smallest such ratio. Similarly, suppose that starting with n copies of the state p'*'^, 
one has a procedure, again involving only local operations and classical communication, for 
distilling m maximally entangled states. The distillable entanglement, D^p"^^), is defined to 
be \0g2d times the asymptotic ratio m/n for an optimal procedure, i.e., one that has the 
largest such ratio. 

A separable state has no entanglement of either sort, whereas a nonseparable state nec- 
essarily has a nonzero entanglement of formation. For pure states, the formation process is 
reversible, so the entanglement of formation and the distillable entanglement are the same. 
For mixed states, however, the distillable entanglement is generally less than the entangle- 
ment of formation, reflecting the irreversibility of the formation process. Interestingly, a 
state with a positive partial transpose has no distillable entanglement [|^. For 2x2 and 
2x3 systems, all entangled states have a nonpositive partial transpose, as noted above, 
and they also have a nonzero distillable entanglement. For 3x3 and higher-dimensional 
systems, however, there are entangled states that have positive partial transpose; though 
these states have nonzero entanglement of formation, they have no distillable entanglement. 
This kind of entanglement, from which no pure-state entanglement can be distilled, is called 
hound entanglement [1^ . 



The entanglement of formation of a pure state \^) of a bipartite system is given by the 
entropy of the marginal density operators, p^ and p'^: 

E{m) = -tr(p^ log2 p^ = -tr(p- log2 p") . (1.3) 

For a bipartite mixed state the entanglement of formation is more complicated. A mixed 
state p'^^ has an ensemble decomposition in terms of pure states l^^j), with probabilities pj, 
if it can be written as 

p"" = Ep.I^.)(^.I- (1-4) 

j 

Modulo a presently unanswered question about the super-additivity of the entanglement 
of formation [^, the entanglement of formation of p"*^ is given by the minimum average 
entanglement of formation of the pure states in an ensemble [0, 



E(p--)= min Vp,E(v[/,), (1.5) 

te.l*j)} j 



where the minimum is taken over all possible ensemble decompositions. For two qubits, 
Wootters has given an explicit formula for the entanglement of formation in terms of what 
he calls the concurrence of the joint density operator p|,p!6| [see Eq. ( |2.4| ) for Wootters's 



concurrence expression]. There are no known explicit formulae for distillable entanglement, 
even for 2x2 systems. We now have all the facts about entanglement we need for posing 
the questions addressed in this paper. 

Let us turn now to issues relevant for distinguishing the theory of real quantum me- 
chanics from standard complex quantum mechanics. The vector space of operators on a 



(i-dimensional vector space (real or complex) is the direct sum of two natural subspaces, the 
space S of real, symmetric matrices, which has dimension ^d{d + 1), and the space A of 
real, antisymmetric matrices, which has dimension ^d{d— 1). If the vector space is over the 
real numbers, all the quantum states and observables lie in the symmetric subspace. For a 
complex vector space, the states and observables are described by Hermitian operators; the 
(real) vector space 7i of Hermitian operators takes advantage of both the symmetric and 
antisymmetric subspaces, it being the direct sum TC = S Q) iA. 

This is of significance for entanglement for the following reason. Suppose one combines 
two systems, with dimensions d^ and dg, to make a composite system with dimension dj^dg. 
In the complex case, the composite vector space of Hermitian operators, of dimension d\d\, 
is the tensor product of the corresponding spaces for A and B, i.e., TIab = 'Ha ® 'Hb- In 
contrast, the composite space of symmetric matrices, of dimension ^d^dB^dAdg + 1), is not 
just the tensor product of the symmetric spaces of A and B, but rather is given by the direct 
sum 

Sab = {Sa ® Sb) © {Aa ® Ab) • (1.6) 

Joint states that have a component in the doubly antisymmetric space, Aa ® Ab, are 
necessarily entangled relative to the real vector space, since product states cannot have a 
component in Aa ® Ab- 

Since any operation in the real vector space can be used in the associated complex vector 
space, an optimal real procedure, either for entanglement of formation E or for distillable 
entanglement D, is never better than an optimal complex procedure. The consequence is 
that for a joint state p^^ in the real vector space, the real (R) and complex (C) entanglement 
measures satisfy the following inequalities: 

Enip-n > Ecipn > Dcipn > DniPn ■ (1-7) 

These ideas are made concrete by considering rebits. The three-dimensional space of 
real, symmetric matrices is spanned by the unit matrix I and two Pauli matrices, a^ and a^, 
whereas the one-dimensional space of real, antisymmetric matrices is spanned by iay. For 
two rebits, the nine- dimensional space Sa ® Sb is spanned by the matrices 

/ ® / cr^0 I az®I 

I® Ox Ox® Ox Oz®Ox , (1.8) 

I ® Oz Ox® Oz Oz® Oz 

but the entire symmetric subspace Sab includes one additional basis matrix, 

iOy ® iOy = —Oy ® Oy . (l.9) 

Any state of the composite system that contains a.Oy®Oy component is necessarily entangled 
relative to the real vector space, simply because tensor products of real states can never sum 
up to give a Oy ® Oy component. 

In Sec. m we show that the entanglement of a two-rebit state is determined entirely by 
the Oy®Oy component of the state: a state p^-^ is separable if and only if tT{p'^'^Oy®Oy) = 0. 
Furthermore, there is a concurrence, defined by 



C(p--) = |tr(p-^V,®a,)|, (1.10) 

which gives the two-rebit entanglement of formation in the same way that Wootters's con- 
currence gives the two-qubit entanglement of formation. In a concluding section (Sec. [m| ), 
we discuss implications of our main result. 

For the present, however, it should be noted that the expression in Eq. ( |1.1U| ) does not 



correspond to Wootters' concurrence formula for qubits simply restricted to real density 
operators. This point is nicely illustrated by the real density operator 

p^s = - (J®/ + ay® ay) . (1.11) 

This state is a separable state relative to complex vector space, as can be checked by the 
partial transpose condition or by noting that it can be derived from an ensemble of two 
product states: 



P"" = \{W + <yy)®W + ^y) + W-<yy)®W- <yy)) 



:i.l2) 



In this (eigen) decomposition, the density operator p"*^ looks like it comes from the mixture 
of two spin states: both particles pointing in the +y direction or both pointing in the —y 
direction. 

In contrast to this, the ensemble decomposition in Eq. ( [L.12| ) is not allowed relative to 
real vector space quantum mechanics. The real concurrence, Eq. ( [L.10| ), of this state is in 



fact C = 1, which means that the state is maximally entangled relative to real vector space. 
Perhaps more interestingly, this state is also a bound entangled state relative to the reals. 
This follows because it is separable relative to the complex numbers, i.e., it has no complex 
entanglement of formation, and hence, by the chain of inequalities in Eq. ( |1.7| ), it has no 
real distillable entanglement. 

II. TWO-REBIT ENTANGLEMENT OF FORMATION 

In this section we first review Wootters's spin-flip operation and how it leads to the 
real concurrence ( |1.10|) , and we then prove our main result, an explicit formula for the real 



entanglement of formation of an arbitrary state of two rebits. 

The spin-flip operation for a single qubit is the anti-unitary operator S = iayC, where C 
denotes complex conjugation in the eigenbasis of az- For a quantum state p of a bipartite 
system — we now drop the superscript AB to reduce clutter in the notation — the spin-flipped 
density operator, distinguished by a tilde, is 

p= {(ry^ay)C{p){cry(^ay) . (2.1) 

The concurrence of a bipartite pure state |\E') is defined to be 

C(vl/)^|(vl/|5|vl/)|. (2.2) 

Defining the concurrence of a mixed state p by 



C{p) = mill Vp,C(^.) , (2.3) 

where the miniinuiii is taken over all possible ensemble decompositions of p, Wootters |^ 
showed that C{p) is given by the explicit expression 



Cw{p) = max [0, ^Ai - ^As - V^3 - ^^4 j , (2.4) 

where the Aj are the (positive) eigenvalues of the operator pp (or of the operator y/ppy/p) 
listed in order of decreasing magnitude. 

For rebit states, the complex conjugation has no effect, and the spin flip simplifies to iCy, 
i.e., a 180° rotation about the y axis. The differences between concurrence and entanglement 
of formation for rebits and qubits can ultimately be traced to the fact that the spin flip for 
qubits is an anti-linear, as opposed to a linear operator. Hence, for a single-rebit pure state 
\ip), the spin- flipped state is 

\^)=tcry\^). (2.5) 

For a joint pure state |\1/) of two rebits, we again define the concurrence to be the overlap 
between |\1/) and the spin-flipped state |\^) = —ay ® ayl"^), i.e. 



and also define the concurrence of a mixed state according to Eq. (2.3). 



The joint pure state |\I^) can be written in terms of a Schmidt decomposition, 

1^) = ailei) ® l/i) + aslea) ® I/2) , (2.7) 

where \ej) and \fj) are the (real) orthonormal eigenvectors of the marginal density operators 
for systems A and B, and ai and 02 are the (positive) square roots of the corresponding 
eigenvalues. It is easy to verify that the concurrence of |\l/) is C(\E') = 2aia2. Thus, as noted 
by Wootters, the concurrence itself can serve as a measurement of entanglement, varying 
smoothly from for product pure states to 1 for maximally entangled pure states. Indeed, 
the entanglement of formation of the pure state |\1/) can be expressed in the form 



Em = H{al) = h(^^^^^^^^ ^ £[C{n , (2.8) 

where 

H{x) = —X log2 X — {1 — x) log2(l — x) (2.9) 

is the binary Shannon entropy. The function S{C) is monotonically increasing and convex 
on the interval < C < 1. 

Before proceeding to the problem at hand, let us delineate a few facts about ensemble 
decompositions of density operators in a real Hilbert space quantum mechanics. Consider a 
mixed state p, defined on a real Hilbert space of dimension d, whose eigendecomposition is 



P = Y.PMi)^^j\ = E |ei)(ej| • (2.10) 

Here the vectors \ej) are the orthonormal eigenvectors of p, with corresponding eigenvalues 
fij, and the vectors \ej) = y/pjlcj) = p^/^|ej) are subnormahzed eigenvectors. In Eq. (p.lO|) 



the sum includes only eigenvectors with nonzero eigenvalues, n < d thus being the rank of 
p. We can restrict attention to the support of p, the subspace of dimension n spanned by 
eigenvectors with nonzero eigenvalue. On this subspace p has a well defined inverse. 
Consider now an arbitrary pure-state ensemble decomposition of p, 

m m 

P = Y.Pj\^i)(^j\ = E ki)K-| • (2-11) 

The ensemble includes m > n normalized vectors \wj), with probabilities pj. The vectors 
\wj) = y/Pj\wj) are subnormahzed, their lengths giving the probabilities. The real version 
of the pure-state decomposition theorem for density operators [p!7|-pi| says that a set of 
subnormahzed vectors gives a decomposition of p if and only if the vectors can be written 

as 



\w 



) = J20kj\ek), j = l,...,m, (2.12) 



k=l 



where O is an n x m matrix whose n rows are real m-dimensional orthonormal vectors. We 
can always extend O to be an m x -m orthogonal matrix by adding additional rows. 
Now notice that the projector onto the support of p can be written as 

m 

^ = Ep-l/2K•)K•|p-l/^ (2.13) 

where 

n 

p-'^'\wj) = Y.Okj\ek), j = l,...,m. (2.14) 

k=l 

By adding additional orthonormal vectors \ek) for /c = n + 1, . . . ,m, we can define m or- 
thonormal vectors in an extended Hilbert space, 

ni 

l^i) - H'^fcjl^fc)' J = l5---,"^- (2.15) 



k=l 



Notice that 



which implies that 



U\wj) = p^^''\wj) , (2.16) 



wj) = p'^'lw,) . (2.17) 



The extension of the resolution of 11 in Eq. ( p.l3|) to a set of orthonormal vectors in 



a 



higher- dimensional space is called a Neumark extension |^ . 



Now we apply these concepts to a two-rebit state p having the pure-state decomposi- 
tion ( p. 11 ). The average concurrence of this decomposition satisfies the inequahty 



(C) = Y.Pfi^'^i) = Y.Pj\i^jWy ® (^y\^j)\ > Y^pA^J 



1(c) I, (2.18) 



where, following Wootters, we define the preconcurrence of the pure state \wj) without the 
absolute value signs that make the concurrence positive: 

c{Wj) = {Wj\(7y ® O'yl'Wj) ■ (2.19) 

The attractive feature of the preconcurrence is that its average value is independent of the 
ensemble decomposition, being given by 

(c) = Y.{^My ® ^y\^j) = E(^J-|P'^'(^J/ ® ('y)p'^^\^j) = tr(r) , (2.20) 

j j 

where 

r = py\ay®ay)p'/'' (2.21) 

is a real, symmetric operator. 

We now show that there is a pure-state ensemble whose average concurrence achieves the 
lower bound, |tr(r)|. We actually show something stronger, using the approach introduced 
by Wootters: there is a pure-state ensemble such that the preconcurrence of every mem- 
ber of the ensemble is tr(r) and, hence, the concurrence of every member of the ensemble 
is |tr(r)|. This stronger result becomes important when we consider the entanglement of 
formation. To construct the desired ensemble, start with the eigendecomposition of p. If 
the eigendecomposition has only one member, then we are dealing with a pure state whose 
preconcurrence is tr(r), and nothing further needs to be done. If the eigendecomposition has 
more than one member, consider the states with the largest and smallest preconcurrences. 
Since the average preconcurrence is tr(r), the largest preconcurrence must be greater than or 
equal to tr(r), and the smallest must be less than or equal to tr(r). Consider the continuous 
sequence of two-dimensional orthogonal matrices that range from the identity matrix to the 
matrix that exchanges the subnormalized vectors for the states with the largest and smallest 
preconcurrences. Somewhere along this sequence, the state with the largest preconcurrence 
is transformed into one whose preconcurrence is tr(r). Adopt the resulting ensemble de- 
composition. Now iterate this procedure, always choosing the states with the largest and 
smallest preconcurrences and transforming the state with the largest preconcurrence to one 
with preconcurrence equal to tr(r). The result is an ensemble decomposition all of whose 
members have preconcurrence equal to tr(r), as promised. 

So far in this section we have shown that for a two-rebit density operator p, the minimum 
average concurrence over all the pure-state ensembles for p is 

{C)^in = niin Y^Pj^i^j) = M^)\ = l^'^iP^y ® ^y)\ ■ (2-22) 



This justifies calling |tr(r)| the concurrence of p, as in Eq. ( 1.10| ) 



The entanglement of formation for a two-rebit state now follows with very little further 
work, since it satisfies the following chain of relations: 

E{p) ^ mxny.PAC{w,)] > min £(j2p,Ciw,)) = £{{0^.) = ^[|tr(r)|] . (2.23) 

The inequality follows from the convexity of the function S{C), and the immediately fol- 
lowing equality follows from the monotonicity of S{C). To saturate the inequality requires 
an ensemble all of whose members have a concurrence equal to (C)min = |tr(r)|. Having 
just constructed an ensemble, we conclude that the entanglement of formation of a two-rebit 
state p is given by 

Eip)=S{MT)\)=S[C{p)]. (2.24) 



III. DISCUSSION 

To conclude, we reiterate that we now possess a complete expression for the entanglement 
of formation of two rebits: 



1 + Jl-C^(p^^ 
E{pn=Hi ^ ^ 1 , (3.1) 

where 

Cipn = Mr)\ = Mp^-ay^ay)\. (3.2) 

In particular, this expression implies that p"*^ is separable relative to real vector space if and 
only if ti {p^'^Cy ^ (Ty) = 0. Notice that this separability condition is equivalent to saying 
that p'^^ is real separable if and only if it is unchanged by partial transposition, i.e., 

pAS ^ ^pAS^T, ^ ^pAS^T, ^ (33) 

where T4 (Tg) denotes partial transposition of system A (B) in any orthonormal basis. 

It is worth stressing the difference between the expression ( p.2|) for the concurrence in 
a real vector space and the Wootters formula (p^.4D for concurrence in standard quantum 



theory. If we let z/j be the eigenvalues of r, ranked in order of decreasing absolute value, the 
real concurrence is 



C(p^^) = |l^l + 1^2 + 1^3 + ^^4 1 . (3.4) 

In contrast, for the Wootters concurrence, one first finds the eigenvalues of tt'^ = r^, these 
being given by A^ = z/|; then the Wootters concurrence is 

Cwip'"'') = max(0, \ui\ - \u2\ - W^l - W4) < C{p^^) . (3.5) 

This difference is illustrated by the class of real states of the form 

9 



pAB 



-{l®I + a{ay®ay)) , (3.6) 



where a is a positive real number that ranges from to 1. For these states, the operator r, 
given by 

T= -(a{I®I) + ay®ay) , (3.7) 

has doubly degenerate eigenvalues \{a ± 1). The real concurrence is C = a, whereas 
the Wootters concurrence is Cw = 0. Thus these states are complex separable, but real 
entangled, except for a = 0. Moreover, because these states are complex separable, they 
have no distillable entanglement, so their real entanglement is bound. 

Our results can be considered a first step toward getting a better understanding of which 
features of quantum entanglement are unique to standard quantum mechanics and which 
are more generic across various foil theories. As just shown, the states (|3.6| ) show that 
bound entanglement is sometimes nothing more than separability with respect to a larger 
field (in this case the complex numbers of standard quantum mechanics). One might ask 
to what extent this is true of bound entanglement in standard quantum mechanics. How 
many bound entangled states in standard quantum mechanics are bound because they are 
separable with respect to a quaternionic theory 0]? 

Another interesting fact is how the regions of entangled vs. separable states within the 
full set of quantum states differ in going from real to complex quantum mechanics. In the 
complex theory, the maximally mixed state p = i/4 of two qubits is surrounded by an open 
set of separable states [^,|^. In the real theory, however, the states (|3.6| ) demonstrate that 
there are entangled states arbitrarily close to the maximally mixed state. 

Where all this will lead, we are not quite sure, but in general it forms part of a larger 
effort to understand the nature of entanglement in our quantum world. 
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